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FOREWORD 


This is a progress report on the research project, “Analysis and Computation of 

Internal Flow Field in a Scramjet Engine,” for the period ended June 30, 1992. Special 
attention during this period was directed to “Investigation of Hypersonic Shock-Induced 
Combustion in a Hydrogen-Air System.” Important results of this study were presented 
at the AIAA 30th Aerospace Sciences Meeting, Reno, Nevada, January 6-9, 1992; AIAA 
paper no. 92-0339, January 1992. This paper is attached in this report as Appendix A. 

This work was supported by the NASA Langley Research Center (Theoretical 
Flow Physics Branch of the Fluid Mechanics Division) through the grant NAG-1-423. 
The grant was monitored by Drs. A. Kumar and J. P. Drummond-Theoretical Flow 
Physics Branch. The work, in part, was also supported by the Old Dominion 
University’s ICAM Program through NASA grant NAG- 1-363; this grant was 
monitored by Mr. Robert L. Yang, Assistant University Affairs Officer, NASA Langley 
Research Center, Hampton, Virginia 23665-5225. 
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INVESTIGATION OF HYPERSONIC SHOCK-INDUCED 


COMBUSTION IN A HYDROGEN-AIR SYSTEM 

J. K. Ahuja* and S. N. Tlwarit 
Old Dominion University, Norfolk, VA 23508 
D. J. Singh* 

Analytical Services and Materials Inc., Hampton, VA 23666 

ABSTRACT 

A numerical study is conducted to simulate the ballistic range experiments at 
Mach 5.11 and 6.46. The flow field is found to be unsteady with periodic instabilities 
originating in the stagnation zone. The unsteadiness of the flow field decreased with 
increase in the Mach number, thus indicating that it is possible to stabilize such flow 
fields with a high degree of overdrive. The frequency of periodic instability is 
determined using Fourier power spectrum and is found to be in good agreement with 
the experimental data. The physics of the instability is explained by the wave interaction 
models available in the literature. 


* Graduate Research Assistant, Department of Mechanical Engineering and Mechanics 
t Eminent Professor, Department of Mechanical Engineering and Mechanics 
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INTRODUCTION 


The national commitment to the National Aerospace Plane (NASP) program and other 
hypersonic vehicles such as Trans-Atmospheric Vehicle (TAV) and Aero-assisted Orbital Transfer 
Vehicle (AOTV) have generated renewed interest in hypersonic flows. Since these vehicles 
will rely on air-breathing propulsion, hypersonic propulsion is one of the key areas actively 
researched. For a successful design of the propulsion system to be used for NASP, it is essential 
to have a clear understanding of the physics of mixing and combustion at supersonic speeds 
in order to develop efficient engines. In the conventional ramjet engine, free-stream air at 
high supersonic speeds is compressed to a low subsonic Mach number at the entrance to the 
combustor. Fuel is injected into the combustor, and burning takes place in a subsonic stream. 
It is advantageous over the standard gas turbines in the Mach number range of 2 to 5, but is 
disadvantageous at hypersonic speeds. Slowing from hypersonic to subsonic speeds will result 
in large pressure losses and will cause very high temperature of air entering the combustor inlet 
(much higher than the adiabatic fuel/air flame temperature), resulting in decomposition of the 
fuel rather than burning. Therefore, the engine will be a drag device rather than a thrust device. 

For an efficient propulsion system at hypersonic speeds, the combustion must take place at 
supersonic speeds, for which two modes of propulsion are being proposed; namely, the Scramjet 
(supersonic combustion ramjet) and Shramjet (shock-induced combustion ramjet). The Scramjet 
([1]— [2]) is an integrated airframe-propulsion concept for a hypersonic airplane. The entire 
undersurface of the vehicle is part of the scramjet engine. Initial compression of the air takes 
place through the bow shock from the nose of the aircraft. Further compression takes place inside 
a series of modules near the rear of the aircraft, thus increasing its pressure and temperature. 
In the combustor, fuel (usually hydrogen) is injected into the hot air by a series of parallel and 
perpendicular injectors where mixing and combustion takes place at supersonic speeds. The 
expansion of burned gases is partially realized through nozzles in the engine modules but mainly 
over the bottom rear surface of the aircraft. At high Mach numbers, the fuel and air do not have 
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enough time for mixing and, therefore, the combustion efficiency decreases. Thus, in order to 
get the desired mixing, the length of the combustor has to be long. Since the highest pressure 
and temperature in the engine occur in the combustor, it has to be very strong; combined with 
the long length, it increases the weight and the drag of the vehicle. 

In order to reduce the size of the combustor, shock- induced combustion (Shramjet [3]) has 
been proposed, where, a shock is employed to increase the temperature of premixed fuel and 
air to a point where chemical reaction will start. Apparent advantages of the Shramjet over 
the Scramjet engine includes very short-length combustors and simple engine geometries. The 
Shramjet’s ability to operate at lower combustor inlet pressures will allow the vehicle to operate 
at a lower dynamic pressure which lessens the heating loads on the airframe. Up to about Mach 
14, Scramjet has better performance than Shramjet, but after that, Shramjet performs better [1]. 

In the past, many researchers have conducted ballistic range experiments to study the 
supersonic combustion/detonation. In these experiments, projectiles were fired in different 
premixed fuel air mixtures, and detonation structures around the projectiles were recorded. Every 
gas mixture has a detonation wave velocity known as Chapman-Jouget (C-J) velocity, which is 
characteristic of the mixture. In any reactive gas mixture, if the normal component of the flow 
velocity following the detonation wave is sonic, then the detonation wave velocity is known 
as the C-J velocity of the mixture. On the other hand, if the normal component of the flow 
velocity is subsonic, the detonation wave velocity is called overdriven, and if supersonic, then it 
is known as underdriven. If the free-stream velocity of the projectile is above the C-J velocity of 
the reactive mixture, the free stream -velocity is referred to as superdetonative. The detonation 
wave structure is highly unstable for projectile velocities less than the Chapman-Jouget velocity 
of the mixture. If the projectile is flying above the C-J velocity of the gas mixture, the detonation 
or reaction front structure shows a coupled shock-deflagration system near the stagnation line of 
the body. These two fronts separate from each other as one moves away from the stagnation 
line. The separation between the two fronts occurs as soon as the velocity component normal to 
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the bow shock is equal to the detonation velocity. The separation between the bow shock and the 
reaction front is called the induction zone. Lehr’s [4] experimental ballistic range shadowgraphs 
for Mach 5.11 and Mach 6.46 are shown in Figs. 1 and 2, respectively. In both cases, a free- 
stream temperature of 292 K and a pressure of 42663.2 N/m 2 (320 mm of Hg) is used along with 
a stoichiometric mixture of hydrogen-air. Fig. 1 shows that the reaction front is separated from 
the bow shock near the stagnation line, and the instabilities or pulsations in the reaction front 
are visible for Mach 5.11. Figure 2 is for Mach 6.46 case, and it is seen that the reaction front 
is coupled with the shock near the stagnation line. A close examination of the shadowgraphs 
reveals that as the flow crosses the bow shock the color changes from light to dark, indicating 
an increase in density. But, as the flow crosses the reaction front, the color changes from dark 
to light, indicating a decrease in density across the reaction front. This is due to a large release 
of energy across the reaction front, causing an increase in the temperature; since the pressure 
remains relatively constant, the density must decrease. 

Behran et al. [5] conducted similar ballistic range experiments by firing plastic spheres into 
hydrogen-air mixtures. Similar behavior of transition from stability to instability of the detonation 
waves was observed when the projectile velocities were decreased to nearly Chapman-Jouget 
velocity of the mixture. 

McVey and Toong [6] also conducted similar experiments where projectiles were fired 
into lean acetylene-oxygen and stoichiometric hydrogen-air mixtures. They proposed the wave 
interaction model to explain the instabilities in the structure of the detonation wave. Their 
model explains how compression waves can be formed when a new reaction front develops in 
the induction zone between the normal segment of the bow shock and the original reaction front. 
These compression waves lead to a cyclic process which is compatible with most of the observed 
features of the flow. However, the strength of the compression waves remained unresolved in 
their wave-interaction model, which is an important factor in determining if such a model is 
physically possible. Alpert and Toong [7] included the effect of the strength of the compression 
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waves and proposed a modified form of the wave-interaction model. 

Several researchers [8-11] have recently attempted to numerically simulate Lehr’s ballistic 
range experiments [4], but have met with limited success. Youngster et al. [8] and Lee and 
Deiwert [9] simulated Lehr’s experimental data for Mach 4.18, 5.11, and 6.46. They used Euler 
equations coupled with species equations to capture the shock and the reaction front. The reaction 
model used was hydrogen-air mixture of six species and an inert gas such as Argon or Nitrogen 
and eight reactions. The flow field was found to be steady in contrast to the experimental 
evidence that the flow field is, indeed, unsteady. For the test conditions of stoichiometric 
hydrogen-air mixture, the detonation wave speed of the mixture is Mach 5.11. Experimentally, 
it has been demonstrated in Lehr’s work that Mach 5.11 and 4.18 show structural instabilities 
of the detonation wave which disappear if the flight Mach number is increased beyond Mach 
5.11. Further, the flow field was not well resolved. They used 32x32 and 57x41 size grids, 
respectively, in their blunt body calculations. These grids were not sufficient to resolve the flow 
field correctly. 

Wilson and MacCormack [10] conducted a detailed numerical investigation of the shock- 
induced combustion phenomena. They used Euler equations and a 13-species and 33-reactions 
chemistry model. They showed the validity of the reaction models and the importance of grid 
resolution needed to properly model the flow physics. They did highly resolved calculations 
for Lehr’s Mach 5.11 and Mach 6.46 cases with adaptive grid. The calculations were not time 
accurate, so that the unsteady behavior was not captured. 

Sussman and Wilson [11] also studied the instabilities in the reaction front for a Mach 
number of 4.79. They also used Euler equations and a 13-species and 33-reactions chemistry 
model. They have proposed a new formulation based on logarithmic transformation. It greatly 
reduces the number of grid points needed to properly resolve the reaction front. They successfully 
simulated the unsteady case. However, the frequency was slightly underpredicted. 

Matsuo and Fujiwara [12] have studied the instabilities of shock-induced combustion around 
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an axisymmetric blunt body. They used Euler equations and a simplified two-step chemistry 
model. They investigated the growth of periodic instabilities by a series of simulations with 
various tip radii and showed that these periodic instabilities are related to shock- standoff distance 
and induction length. They proposed a new model based on McVey and Toong model [6]. The 
instabilities in the reaction front were explained by their model. 

The instability in the structure of the reaction front originates in the induction zone which 
separates the bow shock and the exothermic reaction front in the nose region of the flow field 
and then spreads outwards. In order to capture the physical instabilities, the calculations must 
be carried out for long times to ensure that all relevant time scales are being captured. Since all 
numerical schemes have some numerical diffusion, which is dependent on the grid resolution, a 
coarse grid may damp these oscillations. Further, the numerical damping added to the scheme 
in the vicinity of the reaction front may damp or alter the instability modes. 

The objective of this study is to investigate, in detail, the shock-induced combustion phenom- 
ena for the premixed stoichiometric H 2 -air mixture flow at hypersonic speed (Mach 5.11), which 
is also the Chapman- Jouget speed of the mixture, past a 15 mm spherical projectile, including 
the stability of detonation waves. The analysis is carried out using the axisymmetric version of 
the SPARK2D code [13], which incorporates a seven-species, seven-reactions combustion model 

for hydrogen-air mixture. 
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BASIC GOVERNING EQUATIONS 


The physical model for analyzing the flow field is described by the Navier-Stokes and species 
continuity equations. For two-dimensional axisymmetric flows, these equations are expressed 
in physical coordinates as 


^ + + = u a) 

dt dx dy 


where vectors U, F, G and H are written as 
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The other terms appearing in vectors F, G, and H are defined as 
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In Eq. (1) only (N s -1) species equations need to be considered in the formulation since 
the mass fraction of the species is prescribed by satisfying the constraint equation 

|/, = 1 (12) 
The specific heat at constant pressure for each species is prescribed in Eq. (1 1) by a fourth-order 
polynomial in temperature. The multicomponent diffusion equation for the diffusion velocity 

of the i th species 

U{ = uii + ujj 


is as follows : 
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I, may be noted that this equation has to be applied only to <N,-1) species. The diffusion 
velocity for the remaining species is prescribed by satisfying the constraint equation £ f, u. = 0, 
which ensures the consistency. In Eq.(7>, it is convenient to assume that the body force vector 
per unit mass is negligible. In addition, thermal diffusion is considered to be negligible when 
compared with the binary diffusion coefficient. 
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CHEMISTRY AND THERMODYNAMIC MODELS 

Chemical reaction rate expressions are usually determined by summing the contributions 
from each relevant reaction path to obtain the total rate of change of each species. Each path is 
governed by a law of mass action expression in which the rate constants can be determined from 
a temperature dependent Arrhenius expression. In vector H, the term u>, = M,C, represents the 
net rate of production of species i in all chemical reactions and is modelled as follows . 


N a 


Ns 


^ VjiSi ^ 'ij ~ li ••• 
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1= 1 


u,i = M, Y, Wi - 4 i) 

3 = 1 


N, , N, 
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ro= 1 m = 1 


(15) 


(16) 


where Eq.(14) is a representation of an N r -step chemical reaction, and Eq.(15) is the production 
rate for the i 1 * 1 species as determined from the law of mass action. The reaction constants kq 
and K b j are calculated from the following equations 

■j = ;j = l,. ...N r 07) 
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The equilibrium constant appearing in Eq.(18) is given by 
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( 22 ) 
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The forward rate for each reaction is determined from Eq.(17) which is based on the 
Arrhenius law. The appropriate constants A andej for the H 2 — air reaction system can 
be found in [14], The reverse rate is then calculated from Eq.(18). It should be noted that the 
Gibb’s free energy of each species in g, is obtained from the expression for C p j. 

The hydrogen-air combustion mechanism used in this work is based on the Jachimowski 
Hydrogen-air model [14] which uses seven species and seven reactions. The species are N 2 , O 2 , 
H 2 , OH, H 2 0, O, and H. Each of the seven reactions can proceed in the forward and backward 
directions. The reactions are 

1) 0 2 + H 2 ^OH + OH 

2) 0 2 + H^ OH + O 

3) H 2 + OH ^ H 2 0 + H 

4) H 2 + O ^ OH +H 

5) OH + OH ^ H 2 0 + O 

6) OH + H + N 2 ^ H 2 0 + N 2 

7) H + H + N 2 ^ H 2 + N 2 

The stoichiometric chemical reaction for a hydrogen-air system can be written as 

2H 2 + 0 2 + 3.76N 2 — >2H 2 0 + 3.76N 2 

When a blunt body is moving through a reactive mixture at hypersonic speeds, a bow shock 
is formed ahead of the body, and the temperature of the fuel-air mixture after the bow shock is 
sufficiently high to initiate the reaction. Once the ignition starts, chemical energy is released and 
another discontinuity known as reaction front is formed. In the induction zone, the temperature 
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and the pressure remain relatively constant at the post shock conditions, while the concentrations 
of radicals build up very rapidly. 
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METHOD OF SOLUTION 

The governing equations are transformed from the physical domain (x, y) to a computational 
domain (£, r/) using an algebraic grid generation technique. In the computational domain, Eq. 

(1) is expressed as 


where 


dU_ = 

dt d£ drj 

U = UJ, F = Fy v - Gx v 
G - Gx{ - Fy H = HJ 


(23) 


(24) 


J = x^y v - y^xj, 

The governing equations are solved using the MacCormack’s [17] method. The scheme is 
second-order accurate in time and space. This results in a spatially and temporally discrete, 
simultaneous system of equations at each grid point. The system of equations is solved subject 
to initial and boundary conditions. At the supersonic inflow boundary, all flow quantities are 
specified as free-stream conditions. At the supersonic outflow boundary, all flow quantities are 
extrapolated from interior grid points. Although full N-S equations are used, the slip conditions 
are used to numerically simulate the inviscid flow. A flow tangency or slip boundary condition 
is implied on solid wall. The wall temperature and pressure are extrapolated from interior grid 
points. Initial conditions are obtained by specifying free-stream conditions throughout the flow 
field. The resulting set of equations is marched in time. 

The Lax-Wendroff type schemes are inherently unstable and, hence, higher-order numerical 
dissipation terms are often necessary to get a stable solution. For a non-reacting flow field, an 
artificial viscosity based on temperature and /or pressure is traditionally used, but in chemically 
reacting flows, in addition to temperature and pressure gradients, one can also have very strong 
species concentration gradients. To suppress the numerical oscillations in the induction zone 
where the gradients in the concentration of reactants and products are very strong, additional 
artificial viscosity based on H 2 0 mass fraction is used similar to the one used by Singh et al. [15], 
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RESULTS AND DISCUSSIONS 


The stoichiometric chemical reaction for a hydrogen-air system can be written as 
2H 2 + 0 2 + 3.76N 2 -> 2H 2 0 + 3.76N 2 

When a blunt body is moving through a reactive mixture at hypersonic speeds, a bow shock is 
formed ahead of the body, and the temperature of the fuel-air mixture after the bow shock is 
sufficiently high to initiate the reaction. Once the ignition starts, chemical energy is released and 
another discontinuity known as reaction front is formed. In the induction zone, the temperature 
and the pressure remain relatively constant at the post shock conditions, while the concentrations 
of radicals build up very rapidly. 

Because of the symmetry, only one half of the flow field is calculated. Figure 3 shows 
the typical grid which contains 101x78 grid points (101 normal to the body and 78 along the 
body). For clarity, every fourth grid point is shown in the figure. For the present case of a 
stoichiometric hydrogen-air mixture, the Chapman-Jouget velocity is the same as the velocity of 
the projectile for the Mach 5.11 case. However, for the Mach 6.46 case, the projectile speed 
is significantly above the detonation velocity of the mixture. Unsteady flow phenomenon can 
occur if the free-stream velocity of the projectile is around the C-J detonation velocity of the 
mixture. For both cases, the residuals were dropped by three orders in 12,000 iterations and 
then remained constant. 

Figure 4 shows the contour plot of density for the Mach 5.11 case. The bow shock and 
the reaction front can be seen clearly in the figure. These are separated from each other 
by the induction distance. The separation (i.e., the induction distance) is minimum near the 
stagnation line and increases away from it. This is because near the stagnation line, bow shock 
is almost normal and, hence, the post shock temperature is maximum; thus, induction distance 
is minimum. Away from the stagnation line, the shock strength decreases, thereby decreasing 
the post-shock temperature and, hence, increasing the induction distance. A comparison with 
Fig. 1 shows that all the flow features are very well captured. Figure 5 shows the corresponding 
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plot which has been enlarged for clarity. The bow shock is very crisp and smooth, whereas 
the reaction front is wrinkled. The pulsations or instabilities which arise at the stagnation point 
move through the reaction front. The maximum density is seen to be just after the bow shock, 
and minimum density is after the reaction front. The shock standoff distance is comparable to 
the Lehr’s shadograph. Figure 6 shows the contour plot for temperature, and Fig. 7 shows 
the corresponding enlarged view. Again, it is seen that the bow shock is very smooth, but the 
reaction front which separates from the bow shock near the stagnation line shows pulsations. As 
explained earlier, the post-shock temperature is maximum near the stagnation line and gradually 
decreases away from it. Also, due to the exothermic nature of the reaction, the temperature 
further increases as the reaction proceeds. The peak temperature occurs at the stagnation point. 
Figure 8 shows the temperature along various j=constant grid lines. The post shock stagnation 
point temperature is 3150°K, which compares very well with Ref. [10]. As the gas encounters 
the bow shock, the temperature increases abruptly. Immediately after the shock, the temperature 
stays constant for a short distance and then begins to increase due to exothermic reactions. The 
induction zone decreases with increasing temperature, as chemical energy release will be faster for 
higher temperatures. Also, one can see the unsteadiness in the reaction front. This unsteadiness 
originates from the induction zone near the stagnation line and then travels downstream. The 
contour plots for water mass fraction are shown in Figs. 9 and 10. At the end of the combustion 
zone, the temperature is high enough to start the combustion. As the reaction proceeds, the water 
mass fraction increases rapidly. The oscillations similar to temperature and density profiles can 
be seen here. The instability is characterized by an almost regular periodic wave motion having a 
constant frequency. Similar instability has been observed experimentally in Lehr’s work. Figure 
11 shows the pressure contour. Again, the bow shock is clearly visible in the figure, and the 
pressure jump across the shock is comparable with Ref. [9]. Figure 12 shows the line plot for 
pressure along various j=constant grid lines. As the flow crosses the shock, it encounters the 
pressure jump. The pressure decreases slightly after the shock. The Von Neumann spike, which 
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is characteristic of reacting flows, is also visible. The post shock oscillation in pressure along 
the stagnation line has also been observed in Ref. [9]. 

The instability in the reaction front has been explained by the wave interaction model as pro- 
posed by McVey and Toong [6] and as modified by Matsuo and Fujiwara [11]. Figure 13 shows 
the wave interaction model in terms of the x-t diagram on the stagnation streamline between 
the bow shock and the reaction front. The diagram shows the x-t plot for water mass fraction 
with an overlay of pressure. First, the contact discontinuity approaches the original reaction 
front. The hot gases behind the contact discontinuity begin to react, generating compression 
or pressure waves which propagate upstream and downstream. The compression wave which 
propagates upstream interacts with the bow shock and produces a contact discontinuity behind 
the bow shock. The hot gases on the contact discontinuity behind the bow shock begin to react, 
and thus, generating another set of compression waves. At a somewhat later time, the contact 
discontinuity reaches the position of the original reaction front, extinguishing the reaction at this 
point because no more unreacted gas exists there, and the rate of energy release is effectively 
reduced, and thus, generating rarefaction waves. The reaction front begins to recede because of 
increasing induction time of the colder fluid. The compression wave travelling towards the blunt 
body gets reflected from the body and travels back to the reaction front and causes a change in 
the ignition location, and a new pressure wave is created and then the cycle is repeated. 

Figure 14 shows the x-t plot for density along the stagnation line. The shock front is smooth, 
but the periodic oscillations of the reaction front are clearly visible. These periodic oscillations 
are more clear from Fig. 15, which is the x-t plot for water mass fraction along the stagnation 
line. If one sees these oscillations very closely, it will be clear that the water mass production rate, 
which is also a measure of energy release, continues to increase, and then decreases eventually to 
zero water mass production and, hence, zero energy release. This is the point of extinguishment 
of the reaction front. The reaction almost comes to a standstill at this point. Since the new 
reaction front generated has high energy release (and, hence, high water mass production rate). 
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it sends new sets of compression waves, which propagate both upstream and downstream, and 
the above cycle is repeated. Figure 16 shows the x-t plot for temperature along the stagnation 
line. The periodic oscillations of the reaction front similar to water mass fraction is noticed. 

To further investigate the unsteady nature of the flow field, a Fourier analysis of the flow 
field was conducted. For this, data at various sample stations along the j=61 grid fine were 
stored for 30,000 iterations to get good temporal resolutions. The grid used was 101x78, and 
all calculations were time accurate. Figure 17 shows the amplitude vs frequency plot obtained 
by using Fourier transform. The flow field spectrum is well resolved, and it clearly shows the 
fundamental frequency of 1.2 e+6 Hz and a peak amplitude of 0.004. It also shows subharmonics 
and high-frequency numerical noise. Experimental fundamental frequency, as given in Ref. [16], 
is 1.96e+6 Hz. The discrepancies between the experimental and the numerical value could be due 
to improper grid resolution. The calculations were then repeated for a finer grid (131 x 101). The 
grid aspect ratio was kept the same in both the cases. Figure 18 shows the frequency spectrum 
for the flow field with the finer grid. The sample stations have the same physical locations as 
in the previous case. The dominant frequency now is 2.0e+6 Hz, and the amplitude is 0.004. 
This frequency is in close agreement with the experimental value of 1.96e+6 Hz. The above 
calculations were repeated once again for another finer grid of 197 x 152. The grid aspect ratio 
was kept the same and the sample stations have the same physical locations as in the previous 
cases. Figure 19 shows the frequency spectrum for this grid. The dominant frequency now is 
2.1e+6 Hz., and the amplitude is 0.004. Thus, refining the grid has not changed the frequency 
and therefore, the oscillations in the reaction front are physical. 

The results for the Mach 6.46 case will now be presented. As mentioned earlier, this is a 
superdetonative case, i.e., projectile velocity is higher than the detonation velocity of the mixture. 

The temperature contours are shown in Fig. 20. Qualitatively the results are similar to the 
previous case except near the stagnation zone; the bow shock and the reaction front are now 
almost coupled due to very small induction distance (because of higher post-shock temperature). 
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The two fronts separate from each other slightly downstream. 

Figure 21 shows the water mass-fraction contours. Here, in contrast to the previous case, the 
reaction front is very smooth. The periodic instabilities, which were clearly visible previously, 
cannot be seen. This is in agreement with the experimental result (see Fig. 2). The frequency 
spectrum and the stability for Mach 6.46 has not been analyzed in the present study but will 
be carried on in future work. 
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CONCLUSIONS 


A numerical study is carried out to investigate the shock-induced combustion in premixed 
hydrogen-air mixture. The calculations have been carried out for Mach 5.1 1 and 6.46. The Mach 
5.11 case was found to be unsteady with periodic oscillations. The frequency of oscillations was 
calculated and was found to be in good agreement with the experimentally observed frequency. 
The Mach 6.46 case was found to be macroscopically stable, thus supporting the existing view 
that it is possible to stabilize the shock-induced combustion phenomena with sufficient level 
of overdrive. 
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1‘igmc 2 Shadowgraph of a spherical nose projeclile moving at 
Mach 6.46 into a premixed stoichiometric hydrogen-ail mixture. 
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Figure 5 Enlarged view of density contours for Mach 5. 1 1 



Figure 6 Temperature contours for Mach 5. 1 1 
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Figure 7 Enlarged view of temperature contour for Mach 5. 1 1 
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ABSTRACT 

A numerical study is conducted to simulate the 
ballistic range experiments at Mach 5.11 and 6.46. 
The flow field is found to be unsteady with periodic 
instabilities originating in the stagnation zone. The 
unsteadiness of the flow field decreased with increase 
in Mach number, thus indicating that it is possible 
to stabilize such flow fields with a high degree of 
overdrive. The frequency of periodic instability is 
determined using Fourier power spectrum and is found 
to be in good agreement with the experimental data. 
The physics of the instability are explained by the wave 
interaction model as proposed earlier by McVcy et al. 
and subsequently modified by Matsuo ct al. 

NOMENCLATURE 

Aj reaction rate constant for the j ,h reaction 
Q concentation of i ,!l species 

C p i constant pressure specific heat of i ,h 
species. 

Dj j binary diffusion coefficient of i ,h and j ,h 

species. 

F total (internal and kinetic) energy 
(j activation energy of j 1h reaction 

hj R base enthalpy of i fh species 

k thermal conductivity 

Kfj forward rate constant for j ,h reaction 
Khj backward rate constant for j ,h reaction 
K e qj equilibrium constant for j ,h reaction 

Mj molecular weight of i ,h specie 

N s number of chemical species 

N r number of chemical reactions 

p pressure 

Rj gas constant of i ,h species 

T temperature 

u x-componcnt of the velocity 

t o x-component of the diffusion velocity of 
the i ,h component 

v y-componcnt of the velocity 
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v i y-componcnt of diffusion velocity of 

i ,h species 

Xj mole fraction of i ,h species 
x streamwise coordinate in the physical 
domain 

y normal coordinate in the physical domain 
rvj temperature coefficient in reaction rate 
expression for j ,h reaction 
AGrj Gibbs free energy change for the j th 
reaction 

Anj molar change for the j th reaction 

t] normal coordinate in the computational 
domain 

A second viscosity coefficient 

ft dynamic viscosity 

i/F stoichiometric coefficient of reactant 

corresponding to i ,h species and j lh reaction 
i/' : stoichiometric coefficient of product 

corresponding to i th species and j th reaction 
£ streamwise coordinate in the computational 

domain 
p density 

ct x normal stress in the x-direction 

<7 y normal stress in the y-direction 

r xy shear stress in the xy plane 

u) i production rate of i th species 

INTRODUCTION 

The national commitment to the National Aerospace 
Plane (NASP) program and other hypersonic vehicles 
such as Trans-Almospheric Vehicle (TAV) and Aero- 
assisted Orbital Transfer Vehicle (AOTV) have gener- 
ated renewed interest in hypersonic flows. Since these 
vehicles will rely on air-breathing propulsion, hyper- 
sonic propulsion is one of the key areas actively re- 
searched. For a successful design of the propulsion sys- 
tem to be used for NASP, it is essential to have a clear 
understanding of the physics of mixing and combus- 
tion at supersonic speeds in order to develop efficient 
engines. In the conventional ramjet engine, free-stream 
air at high supersonic speeds is compressed to a low 
subsonic Mach number at the entrance to the combus- 
tor. Fuel is injected into the combustor, and burning 
takes place in a subsonic stream. It is advantageous 
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over the standard gas turbines in the Mach number 
range of 2 to 5, but is disadvantageous at hypersonic 
speeds. Slowing from hypersonic to subsonic speeds 
will result in large pressure losses and will cause very 
high temperature of air entering the combustor inlet 
(much higher than the adiabatic fuel/air flame tempera- 
ture), resulting in decomposition of the fuel rather than 
burning. Therefore, the engine will be a drag device 
rather than a thrust device. 

For an efficient propulsion system at hypersonic 
speeds, the combustion must take place at supersonic 
speeds, for which two modes of propulsion are be- 
ing proposed; namely, the Scramjet (supersonic com- 
bustion ramjet) and Shramjet (shock-induced combus- 
tion ramjet). The Scramjet ([ 1 ]-f2J) is an integrated 
airframe-propulsion concept for a hypersonic airplane. 
The entire undersurface of the vehicle is part of the 
scramjet engine. Initial compression of the air takes 
place through the bow shock from the nose of the air- 
craft. Further compression takes place inside a series of 
modules near the rear of the aircraft, thus increasing its 
pressure and temperature. In the combustor, fuel (usu- 
ally hydrogen) is injected into the hot air by a series 
of parallel and perpendicular injectors where mixing 
and combustion takes place at supersonic speeds. The 
expansion of burned gases is partially realized through 
nozzles in the engine modules but mainly over the bot- 
tom rear surface of the aircraft. At high Mach num- 
bers, the fuel and air do not have enough time for 
mixing and, therefore, the combustion efficiency de- 
creases. Thus, in order to get the desired mixing, the 
length of the combustor has to be long. Since the high- 
est pressure and temperature in the engine occur in the 
combustor, it has to be very strong; combined with the 
long length, it increases the weight and the drag of the 
vehicle. 

In order to reduce the size of the combustor, 
shock-induced combustion (Shramjet [31) has been pro- 
posed, where, a shock is employed to increase the tem- 
perature of premixed fuel and air to a point where 
chemical reaction will start. Apparent advantages of 
the Shramjet over the Scramjet engine includes very 
short-length combustors and simple engine geometries. 
The Shramjet’ s ability to operate at lower combustor in- 
let pressures will allow the vehicle to operate at a lower 
dynamic pressure which lessens the heating loads on 
the airframe. Up to about Mach 14, Scramjet has bet- 
ter performance than Shramjet, but after that, Shramjet 
performs better [1]. 

In the past, many researchers have conducted bal- 
listic range experiments to study the supersonic com- 
bustion/detonation. In these experiments, projectiles 
were fired in different premixed fuel air mixtures, 
and detonation structures around the projectiles were 
recorded. Every gas mixture has a detonation wave ve- 
locity known as Chapman-Jouget (C-J) velocity, which 


is characteristic of the mixture. In any reactive gas 
mixture, if the normal component of the flow velocity 
following the detonation wave is sonic, then the deto- 
nation wave velocity is known as the C-J velocity of 
the mixture. On the other hand, if the normal com- 
ponent of the flow velocity is subsonic, the detonation 
wave velocity is called overdriven, and if supersonic, 
then it is known as underdriven. If the free-stream ve- 
locity of the projectile is above the C-J velocity of the 
reactive mixture, the free stream-velocity is referred to 
as supcrdetonalive. The detonation wave structure is 
highly unstable for projectile velocities less titan the 
Chapman-Jouget velocity of the mixture. If the projec- 
tile is flying above the C-J velocity of the gas mixture, 
the detonation or reaction front structure shows a cou- 
pled shock-deflagration system near the stagnation line 
of the body. These two fronts separate from each other 
as one moves away from the stagnation line. The sep- 
aration between the two fronts occurs as soon as the 
velocity component normal to the bow shock is equal 
to the detonation velocity. The separation between die 
bow shock and the reaction front is called the induction 
zone. Lehr's [4] experimental ballistic range shadow- 
graphs for Mach 5. 1 1 and Mach 6.46 are shown in Figs. 

1 and 2, respectively. In both cases, a free-stream tem- 
perature of 292 K and a pressure of 42663.2 N/m 2 (320 
mm of Hg) is used along with a stoichiometric mixture 
of hydrogen-air. Fig. 1 shows that the reaction front is 
separated from the bow shock near the stagnation line, 
and the instabilities or pulsations in the reaction front 
are visible for Mach 5.11. Figure 2 is for Mach 6.46 
case, and it is seen that the reaction front is coupled 
with the shock near the stagnation line. A close ex- 
amination of the shadowgraphs reveals that as the flow 
crosses the bow shock the color changes from light 
to dark, indicating an increase in density. But, as the 
flow crosses the reaction front, the color changes from 
dark to light, indicating a decrease in density across 
the reaction front. This is due to a large release of 
energy across the reaction front, causing an increase in 
the temperature, since the pressure remains relatively 
constant, the density must decrease. 

Behran et al. [5) conducted similar ballistic range 
experiments by firing plastic spheres into hydrogen-air 
mixtures. Similar behavior of transition from stability 
to instability of the detonation waves was observed 
when the projectile velocities were decreased to nearly 
Chapman-Jouget velocity of the mixture. 

McVey and Toong (6| also conducted similar 
experiments where projectiles were fired into lean 
acetylene-oxygen and stoichiometric hydrogen-air mix- 
tures. They proposed the wave interaction model to 
explain the instabilities in die structure of the deto- 
nation wave. Their model explains how compression 
waves can be formed when a new reaction front devel- 
ops in the induction zone between the normal segment 
of the bow shock and the original reaction front. These 
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compression waves lead to a cyclic process which is 
compatible with most of the observed features of the 
How. However, the strength of the compression waves 
remained unresolved in their wave-interaction model, 
which is an important factor in determining if such a 
model is physically possible. Alperl and Toong [7] in- 
cluded the effect of the strength of the compression 
waves and proposed a modified form of the wave- 
interaction model. 

Several researchers [8-1 1 ] have recently attempted 
to numerically simulate Lehr's ballistic range experi- 
ments [4], but have met with limited success. Young- 
ster et al. [8] and Lee ct al. [9| simulated Lehr’s 
experimental data for Mach 4.18, 5.11, and 6.46. Ihcy 
used Huler equations coupled with species equations to 
capture the shock and the reaction front. The reaction 
model used was hydrogen-air mixture of six species 
and an inert gas such as Argon or Nitrogen and eight 
reactions. The flow field was found to be steady in con- 
trast to the experimental evidence that the flow field 
is, indeed, unsteady. For the test conditions of stoi- 
chiometric hydrogen-air mixture, the detonation wave 
speed of the mixture is Mach 5.11. Hxperimentally, it 
has been demonstrated in Lehr's work that Mach 5.11 
and 4.18 show structural instabilities of the detonation 
wave which disappear if the (light Mach number is 
increased beyond Mach 5.11. Further, the flow field 
was not well resolved. Ihey used 32x32 and 57x41 
size grids, respectively, in their blunt body calculations. 
These grids were not sufficient to resolve the How field 
correctly. 

Wilson ct al. [101 conducted a detailed numeri- 
cal investigation of the shock-induced combustion phe- 
nomena. They used Euler equations and a 13-species 
and 33-rcactions chemistry model, lliey showed the 
validity of the reaction models and the importance 
of grid resolution needed to properly model the (low 
physics, lliey did highly resolved calculations for 
Lehr's Mach 5.11 and Mach 6.46 cases with adaptive 
grid. The calculations were not time accurate, so that 
the unsteady behavior was not captured. 

Sussman et al. [11) also studied the instabilities 
in the reaction front for a Mach number of 4.79. They 
also used fuller equations and a 13-spccics and 33-re- 
actions chemistry model. They have proposed a new 
formulation based on logarithmic transformation. It 
greatly reduces the number of grid points needed to 
properly resolve the reaction front, lliey successfully 
simulated the unsteady case. However, the frequency 
was slightly undcrpredicted. 

Matsuo and Eujiwara [12| have studied the in- 
stabilities of shock-induced combustion around an ax- 
isymmctric blunt body. lliey used Euler equations and 
a simplified two-step chemistry model. They investi- 
gated the growth of periodic instabilities by a scries 
of simulations with various tip radii and showed that 


these periodic instabilities are related to shock-standoff 
distance and induction length. They proposed a new 
model based on McVcy and Toong's model [61. The in- 
stabilities in the reaction front were explained by their 
model. 


The instability in the structure of the reaction front 
originates in the induction zone which separates the 
bow shock and the exothermic reaction front in the 
nose region of the flow field and then spreads outwards. 
In order to capture the physical instabilities, the calcu- 
lations must be carried out for long times to ensure 
that all relevant time scales are being captured. Since 
all numerical schemes have some numerical diffusion, 
which is dependent on the grid resolution, a coarse 
grid may damp these oscillations. Further, the numeri- 
cal damping added to the scheme in the vicinity of the 
reaction front may damp or alter the instability modes. 
The objective or this study is to investigate, in de- 
tail, the shock-induced combustion phenomena for the 
premixed stoichiometric l^-air mixture flow at hyper- 
sonic speed (Mach 5.11), which is also the Chapman- 
Jouget speed of the mixture, past a 15 mm spherical 
projectile, including the stability of detonation waves. 
The analysis is carried out using the axisymmetric ver- 
sion of the SPARK2D code [13], which incorporates 
a seven-species, seven-reactions combustion model for 
hydrogen-air mixture. 

BASIC GOVERNING EQUATIONS 

The physical model for analyzing the flow field is 
described by the Navier-Stokes and species continuity 
equations. For two-dimensional axisymmctric flows, 
these equations are expressed in physical coordinates 
as 


dV dF_ (]G_ 

Ot dr dy 
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where vectors II, F, G and II are written as 
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In Eq. (l)only (N s — 1) species equations need to 
be considered in the formulation since the mass fraction 
of the species is prescribed by satisfying the constraint 
equation 

n. 

£/i = 1 02) 

i = 1 

The specific heat at constant pressure for each species 
is prescribed in Hq. (1 1) by a fourth-order polynomial 
in temperature. The multicomponent diffusion equation 
for the diffusion velocity of the i ,h species 

u, = Uii + tijj (13) 


is as follows : 


The other terms appearing in vectors F f G, and II 
are detined as 
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It may be noted that this equation has to be applied 
only to (N* — 1) species, lhe diffusion velocity for the 

remaining species is prescribed by satisfying the con* 
N. 

slraint equation h — 0, which ensures the con- 

1=1 

sistency. In Eq.(7), it is convenient to assume that the 
body force vector per unit mass is negligible. In ad- 
dition, thermal diffusion is considered to be negligible 
when compared with the binary diffusion coefficient. 

CHEMISTRY AND 
THERMODYNAMIC MODELS 

Chemical reaction rate expressions are usually de- 
termined by summing the contributions from each rel- 
evant reaction path to obtain the total rate of change 
of each species. Each path is governed by a law of 
mass action expression in which the rale constants can 
be determined from a temperature dependent Arrhe- 
nius expression. In vector II, the term u ) t = M t C\ 
represents the net rate of production of species i in all 
chemical reactions and is modelled as follows : 


N. 

N, 





1 = 1 

11 

EM 

Si ; j = 1, 

...N r 

(15) 

N r 


r 

N. 

N. 

a 

£ K - 

- *U) 

K I , 

[j rtr 

n 

C , '>« 

1 r/i 

j- \ 


_ 

rn ~ 1 

m = 1 







(16) 


where Eq.(14) is a representation of an N r -step chemi- 
cal reaction, and Eq.(15) is the production rate for the 
i ,h species as determined from the law of mass action. 
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The reaction constants Kfj and K b j are calculated from 
the following equations 


rj) using an algebraic grid generation technique. In the 
computational domain, Hq. (1) is expressed as 


«„ = A,T*«*(0) :i = > Hr C7) 


au_ ar 

Ot + di] 


= H 


(23) 


Kbj 


i = = 1 Nr 


(18) 




The equilibrium constant appearing in Iiq.( 18) is 
given by 



where 


N, 


i=l 


A”; =!>"<-£ *4 ‘ = 

i = l 


= 1, ... A^ r (20) 




N. 


* (21) 


iL = A,-T(l-/nT)- 
Ri 




(tI) 1 " - (I) 1 * + 


Fi-CiT;i= I Af, 

( 22 ) 


The forward rate for each reaction is determined 
from Hq.(17) which is based on the Arrhenius law. 1 he 
appropriate constants A j,«j, and for the Il 2 — air 
reaction system can be found in (14). The reverse rate 
is then calculated from F!q.(18). It should be noted that 
the Gibb's free energy of each species in gi is obtained 
from the expression for C pi . 

The hydrogen-air combustion mechanism used in 
Ibis work is based on the Jachimowski Hydrogen-air 
model [14] which uses seven species and seven reac- 
tions. The species are N 2 , O 2 . H 2 , H 2 D, D, and II. 
I?ach of the seven reactions can proceed in the forward 
and backward directions. The reactions are 

1) () 2 + H 2 —OH + OH 

2) 0 2 + II— OH + O 

3) 1I 2 + OH — H 2 0 + II 

4) 1I 2 + O =! OH +H 

5) Oil + Oil — II 2 O + o 

6) OH + II + N 2 — IhO + N 2 

7) II + 11 + N 2 — H 2 + N 2 

METHOD OF SOLUTION 

The governing equations are transformed lrom the 
physical domain (x, y) to a computational domain (£, 


where 

U = UJ, F= Fy n -Gx n 
G = Gx ( -Fyt,H = HJ < 24 > 

./ = XtV n - y^r, 

The governing equations are solved using MacCor- 
mack's [171 method. Ilic scheme is second-order ac- 
curate in time and space. This results in a spatially 
and temporally discrete, simultaneous system of equa- 
tions at each grid point. The system of equations is 
solved subject to initial and boundary conditions. At 
the supersonic inflow boundary, all flow quantities arc 
specified as frcc-stream conditions. At the supersonic 
outflow boundary, all flow quantities arc extrapolated 
from interior grid points. Although full N-S equations 
are used, the slip conditions are used to numerically 
simulate the inviscid flow. A flow tangency or slip 
boundary condition is implied on solid wall. Hie wall 
temperature and pressure are extrapolated from interior 
grid points. Initial conditions are obtained by speci- 
fying free-slream conditions throughout the flow field. 
The resulting set of equations is marched in time. 

llte I.ax-Wcndroff type schemes arc inherently un- 
stable and, hence, higher-order numerical dissipation 
terms arc often necessary to get a stable solution. I-or 
a non-reacting flow field, an artificial viscosity based 
on temperature and /or pressure is traditionally used, 
but in chemically reacting flows, in addition to temper- 
ature and pressure gradients, one can also have very 
strong species concentration gradients. Io suppress the 
numerical oscillations in the induction zone where the 
gradients in the concentration of reactants and products 
are very strong, additional artificial viscosity based on 
I1 2 0 mass fraction is used similar to the one used by 
Singh et al. [151 

RESULTS AND DISCUSSIONS 

The stoichiometric chemical reaction for a 
hydrogen-air system can be written as 

2H 2 + 0 2 + 3.76N 2 — 2H 2 0 + 3.76N 2 
When a blunt body is moving through a reactive mix- 
ture at hypersonic speeds, a bow shock is formed ahead 
of the body, and the temperature of the fuel-air mixture 
after the bow shock is sufficiently high to initiate the 
reaction. Once the ignition starts, chemical energy is 
released and another discontinuity known as reaction 
front is formed. I 11 the induction zone, the temperature 
and the pressure remain relatively constant at the post 
shock conditions, while the concentrations of radicals 
build up very rapidly. 
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Because of the symmetry, only one half of the 
flow field is calculated. Figure 3 shows the typical 
grid which contains 101x78 grid points (101 normal 
to the body and 78 along the body). For clarity, ev- 
ery fourth grid point is shown in the figure. For the 
present case of a stoichiometric hydrogen-air mixture, 
the Chapman-Jouget velocity is the same as the velocity 
of the projectile for the Mach 5.11 case. However, for 
the Mach 6.46 case, the projectile speed is significantly 
above the detonation velocity of the mixture. Unsteady 
flow phenomenon can occur if the free-stream velocity 
of the projectile is around the C-J detonation velocity of 
the mixture. For both cases, the residuals were dropped 
by three orders in 12,000 iterations and then remained 
constant. Figure 4 shows the contour plot of density 
for the Mach 5.11 case. The bow shock and the reac- 
tion front can be seen clearly in the figure. They are 
separated from each other by the induction distance. 
The separation (i.e., the induction distance) is mini- 
mum near the stagnation line and increased away from 
it. This is because near the stagnation line, bow shock 
is almost normal and, hence, the post shock tempera- 
ture is maximum; thus, induction distance is minimum. 
Away from the stagnation line, the shock strength de- 
creases, thereby decreasing the post-shock temperature 
and, hence, increasing the induction distance. A com- 
parison with Fig. 1 shows that all the flow features are 
very well captured. Figure 5 shows the correspond- 
ing plot which has been enlarged for clarity. The bow 
shock is very crisp and smooth, whereas the reaction 
front is wrinkled. The pulsations or instabilities which 
arise at the stagnation point move all through the re- 
action front. The maximum density is seen to be just 
after the bow shock, and minimum density is after the 
reaction front. The shock standoff distance is compa- 
rable to the Lehr’s shadograph. Figure 6 shows the 
contour plot for temperature, and Fig. 7 shows the 
corresponding enlarged view. Again, it is seen that 
the bow shock is very smooth, but the reaction front 
which separates from the bow shock near the stagna- 
tion line shows pulsations. As explained earlier, the 
post-shock temperature is maximum near the stagna- 
tion line and gradually decreases away from it. Also, 
due to the exothermic nature of the reaction, the tem- 
perature further increases as the reaction proceeds. The 
peak temperature occurs at the stagnation point. Fig- 
ure 8 shows the temperature along various j=constant 
grid lines. The post shock stagnation point temperature 
is 3 1 50° K, which compares very well with Ref. [10]. 
As the gas encounters the bow shock, the temperature 
increases abruptly. Immediately after the shock, the 
temperature stays constant for a short distance and then 
begins to increase due to exothermic reactions. The in- 
duction zone decreases with increasing temperature, as 
chemical energy release will be faster for higher tem- 
peratures. Also, one can see the unsteadiness in the 
reaction front. This unsteadiness originates from the 


induction zone near the stagnation line and then travels 
downstream. The contour plots for water mass fraction 
are shown in Figs. 9 and 10. At the end of the combus- 
tion zone, the temperature is high enough to start the 
combustion. As the reaction proceeds, the water mass 
fraction increases rapidly. The oscillations similar to 
temperature and density profiles can be seen here. The 
instability is characterized by an almost regular peri- 
odic wave motion having a constant frequency. Similar 
instability has been observed experimentally in Lehr’s 
work. Figure 1 1 shows the pressure contour. Again, 
the bow shock is clearly visible in the figure, and the 
pressure jump across the shock is comparable with Ret. 
[9]. Figure 12 shows the line plot for pressure along 
various j=constant grid lines. As the flow crosses the 
shock, it encounters the pressure jump. The pressure 
decreases slightly after the shock, llie Von Neumann 
spike, which is characteristic of reacting flows, is also 
visible. 'Hie post shock oscillation in pressure along 
the stagnation line has also been observed in Ret |9|. 

The instability in the reaction front has been ex- 
plained by the wave interaction model as proposed by 
McVey and Toong [6] and as modified by Matsuo and 
Fujiwara [11]. Figure 13 shows the wave interaction 
model in terms of the x-t diagram on the stagnation 
streamline between the bow shock and the reaction 
front. The diagram shows the x-l plot for water mass 
fraction with an overlay of pressure. First, the con- 
tact discontinuity approaches the original reaction front. 
'Hie hot gases behind the contact discontinuity begin to 
react, generating compression or pressure waves which 
propagate upstream and downstream. Ihe compression 
wave which propagates upstream interacts with the bow 
shock and produces a contact discontinuity behind the 
bow shock. The hot gases on the contact discontinuity 
behind the bow shock begin to react, and thus, generat- 
ing another set of compression waves. At a somewhat 
later lime, the contact discontinuity reaches the posi- 
tion of the original reaction front, extinguishing the 
reaction at this point because no more unreacted gas 
exists there, and the rate of energy release is effec- 
tively reduced, and thus, generating rarefaction waves. 
The reaction front begins to recede because of increas- 
ing induction time of the colder fluid. The compression 
wave travelling towards the blunt body gets reflected 
from the body and travels back to the reaction trout 
and causes a change in the ignition location, and a new 
pressure wave is created and then the cycle is repeated. 

Figure 14 shows the x-t plot for density along the 
stagnation line. The shock front is smooth, but the peri- 
odic oscillations of the reaction front are clearly visible. 
These periodic oscillations are much more clear from 
Fig. 15, which is the x-t plot for water mass fraction 
along the stagnation line. If one sees these oscillations 
very closely, it will be clear that the water mass pro- 
duction rate, which is also a measure of energy release, 
continues to increase, and then decreases eventually to 
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zero water mass production and, hence, zero energy 
release. This is the point of extinguishment of the re- 
action front. The reaction almost comes to a standstill 
at this point. Since the new reaction front generated has 
high energy release (and, hence, high water mass pro- 
duction rate), it sends new sets of compression waves, 
which propagate hoth upstream and downstream, and 
the above cycle is repeated. Figure 16 shows the x-t 
plot for temperature along the stagnation line. The pe- 
riodic oscillations of the reaction front similar to water 
mass fraction is noticed. 

To further investigate the unsteady nature of the 
llow field, a Fourier analysis of the (low field was con- 
ducted. For this, data at various sample stations along 
the j=61 grid line were stored for 30,000 iterations 
to get good temporal resolutions. The grid used was 
101x78, and all calculations were time accurate. Fig- 
ure 17 shows the amplitude vs frequency plot obtained 
by using Fourier transform. The flow field spectrum 
is well resolved, and it clearly shows the fundamen- 
tal frequency of 1.2 e+6 Hz and a peak amplitude of 

0.004. It also shows subharmonics and high-frequency 
numerical noise. Experimental fundamental frequency, 
as given in Ref. (16], is 1.96e+6 Hz. The discrepan- 
cies between the experimental and the numerical value 
could be due to improper grid resolution. The calcu- 
lations were then repeated for a finer grid (131x101). 
Ilic grid aspect ratio was kept the same in both the 
cases. Figure 18 shows the frequency spectrum for the 
fiow field with the liner grid. The sample stations have 
the same physical locations as in the previous case. 'Hie 
dominant frequency now is 2.0e+6 Hz, and the ampli- 
tude is 0.004. This frequency is in close agreement 
with the experimental value of l.96e+6 Hz. The above 
calculations were repeated once again for another finer 
grid of 197x152. The grid aspect ratio was kept the 
same and the sample stations have the same physical 
locations as in the previous cases. Figure 19 shows the 
frequency spectrum for this grid Ihc dominant fre- 
quency now is 2. 1 e+6 Hz., and the amplitude is 0.004. 
Hms, refining the grid has not changed the frequency 
and therefore, the oscillations in the reaction front arc 
physical. 

Hie results lor the Mach 6.46 case will now be 
presented. As mentioned earlier, this is a superdeto- 
native case, i.e„ projectile velocity is higher than the 
detonation velocity of the mixture 

The temperature contours are shown in Fig. 20. 
Qualitatively the results are similar to the previous case 
except near the stagnation zone; the bow shock and 
the reaction front arc now almost coupled due to very 
small induction distance (because of higher post-shock 
temperature). The two fronts separate from each other 
slightly downstream. 

Figure 21 shows the water mass-fraction contours. 
Here, in contrast to the previous case, the reaction front 


is very smooth. The periodic instabilities, which were 
clearly visible previously, cannot be seen This is in 
agreement with the experimental result (see Fig. 2). 
'Hie frequency spectrum and the stability for Mach 6.46 
has not been analyzed in the present study but will be 
carried on in future work. 

CONCLUSIONS 

A numerical study is carried out to investigate 
the shock-induced combustion in premixed hydrogen- 
air mixture. The calculations have been carried out for 
Mach 5.1 1 and 6.46. The Mach 5.1 1 case was found to 
be unsteady with periodic oscillations. 'Hie frequency 
of oscillations was calculated and was found to be in 
good agreement with the experimentally observed fre- 
quency. The Mach 6.46 case was found to be macro- 
scopically stable, thus supporting the existing view that 
it is possible to stabilize the shock-induced combustion 
phenomena with sufficient level of overdrive. 
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Figure I Shadowgraph of a spherical nose projectile moving at 
Mach 5.1 I into a premixed stoichiometric hydrogen-air mixture. 



Figure 2 Shadowgraph of a spherical nose projectile moving at 
Mach 6.46 into a premixed stoichiometric hydrogen-air mixture. 
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Figure 3 Typical grid for the calculation 
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Figure 5 Enlarged view of density 
contours for Mach 5.11 




Figure 4 Density contours for Mach 5. 1 1 Figure 6 Temperature contours for Mach 5. 1 1 
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Figure 7 Enlarged view of temperature 
contour for Mach 5.1 1 



Figure 8 Temperature vs radius along 
various j= constant grid lines 



Figure 9 Water mass fraction 
contours for Mach 5.11 
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Figure 10 Enlarged view of H 2 O 
mass fraction contours 
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Figure 1 1 Pressure contours tor Mach 5.1 1 



Figure 13 x-t plot for the wave 
interaction model (H 2 0 contours 
with an overlay of pressure) 




Figure 12 Pressure vs radius along Figure 14 x-t plot tor density 

various j= constant grid lines along the stagnation line 
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Figure 15 x-t plot for H 2 O mass-fraction 
along stagnation line 



Figure 17 Temporal frequency spectrum of 
water mass fraction for 101x78 grid size. 



x 

Figure 16 x-t plot for temperature 
along stagnation line 



Figure 18 Temporal frequency spectrum of 
water mass fraction for 131x101 grid size. 
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